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We study the quantum dynamics generated by the two-axis counter-twisting Hamiltonian from an 
initial spin coherent state in a spin-1/2 ensemble. A characteristic feature of the two-axis counter¬ 
twisting Hamiltonian is the existence of four neutrally stable and two saddle unstable fixed points. 
The presence of the last one is responsible for a high level of squeezing. The squeezing is accompanied 
by the appearance of several quantum states of interest in quantum metrology with Heisenberg- 
limited sensitivity, and we show fidelity functions for some of them. We present exact results for 
the quantum Fisher information and the squeezing parameter. Although, the overall time evolution 
of both changes strongly with the number of particles, we find that they have regular dynamics for 
short times. We explain scaling with the system size by using a Gaussian approach. 

PACS numbers: 03.75.Gg, 03.75.Dg. 


I. INTRODUCTION 

Kitagawa and Ueda in their pioneering work pQ have 
proposed the one-axis twisting (OAT) and two-axis 
counter-twisting (TACT) Hamiltonians for dynamical 
generation of spin-squeezed states. The Bose-Einstein 
condensate of ultra-cold atoms offers an exceptional 
tool for experimental realization of the one-axis twisting 
Hamiltonian and investigation of the many-particle en¬ 
tanglement. Several experiments have reported the cre¬ 
ation of spin-squeezed states, by manipulation of inter¬ 
nal states of multicomponent condensates US, or alter¬ 
natively in a single-component condensate in a double¬ 
well potential El IS. In contrary, the two-axis counter¬ 
twisting Hamiltonian cannot be simply realized by inter¬ 
atomic interactions among ultra-cold atoms. Therefore 
several schemes were proposed to transform the one-axis 
twisting Hamiltonian into an effective two-axis counter¬ 
twisting Hamiltonian HHn] or implement TACT model 
in other realistic systems CM!]. Actually, there is grow¬ 
ing interest in quantum states generated by the TACT 
Hamiltonian since they support Heisenberg-like sensitiv¬ 
ity of high-precision measurements and a much higher 
level of squeezing that is unachievable by OAT interac¬ 
tions EG!!- 

In this paper, we study in great detail the quan¬ 
tum dynamics generated by the two-axis counter-twisting 
Hamiltonian from an initial spin coherent state in a spin- 
1/2 ensemble. We start with the mean-field description 
which identifies a convenient location of the initial coher¬ 
ent spin state. A characteristic feature of the two-axis 
counter-twisting Hamiltonian is the existence of the four 
neutrally stable and two saddle unstable fixed points. 
The presence of the last one is responsible for high level 
of squeezing. On the other hand, the quantum dynamics 
around a stable fixed point generates states that support 
shot-noise limited sensitivity for quantum metrology. 

The location of the initial spin coherent state on an un¬ 
stable fixed point leads to strong stretching of the state 
along a meridian resulting in highly reduced variance of 
the spin operator. It has to be noted that the angle 


between the inflowing and outflowing trajectories in the 
phase space is 7r/2, thus the squeezing is optimal. The 
high level of squeezing is accompanied by the presence of 
several quantum states of interest in quantum metrology 
with Heisenberg-limited sensitivity. We will show high 
fidelity to the Berry-Wiseman, equally-weighted super¬ 
position and Yurke states, and a bit worse to the twin- 
Fock state, as it was partially reported in mi- We show 
that better sensitivity can be reached by using states pro¬ 
duced by the dynamics at later times than the time for 
the optimal squeezing. 

We calculate the quantum Fisher information in order 
to quantify the amount of quantum correlations gener¬ 
ated in time that are useful for precision measurements. 
Although, the overall time evolution of the quantum 
Fisher information changes strongly when one changes 
the number of particles N, we find that it has regular 
dynamics for short times of interest. In this regime, the 
time scales like ~ ln(27rA) /N, and we explain the scaling 
using a Gaussian approach within the Bogoliubov-Born- 
Green-Kirkwood-Yvon hierarchy mm, as well as the 
known scaling of the variance of the spin operator. Our 
results show that the quantum dynamics with the two- 
axis counter-twisting Hamiltonian creates quantum cor¬ 
relations in a regular way on a short time scale. This time 
scale is reduced with an increased number of particles. 


II. THE MODEL 

We consider a collection of N qubits e.q. particles in 
two modes. The system is conveniently described using 

the collective spin operator S , whose components written 
in the Schwinger representation are 

S x = \ + b f a^j , (la) 

S y = 2 - tfa) , (lb) 

S z = 7 } - tfb'j , (lc) 
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where a), w are creation operators associated with two 
modes. The two-axis counter-twisting Hamiltonian pro¬ 
posed by Kitagawa and Ueda |T| is 

^TACT = §(S+- S-) . (2) 

where S± = S x ± iS'y. Nevertheless, we will operate 
on the rotated Hamiltonian H = with U = 

e -z(7T/2)S B) 

n = -hx ( s v s z + s z s y ). (3) 

Although, the action of a SU(2) group element on <© 
does not change the overall result, it describes a distinct 
physical system. In our case, the Hamiltonian ^ sim¬ 
plifies the form of observables of interest. 

The Schrodinger equation 

mMt)) = u\*(t)) (4) 


cannot be solved analytically. We solve it numerically in 
the Fock state basis with fixed number of particles N. 
A pure state can be decomposed in the Fock state basis 

N 

|\Ii(t)) = Y c k{t) | k, N — k), and then coupled first-order 
k =o 

differential equations for the coefficients Ck can be solved 
using matrix exponential. In this way, one can forget 
about an abstract Hilbert space and consider only the 
representation of operators and states in a more familiar 
vector space. Operations such as dot product, addition 
and multiplication transfer into the vector space. Ob¬ 
servables are represented as square hermitian matrices, 
and ket states as column vectors. 

The initial state for the evolution is a spin coherent 
state written in the Fock basis (18| . 


N 

\ e ^) = 


k =0 




x | k,N — k), (5) 

and is parametrized by the two real variables 0 < 9 < tt 
and —tt < ip < tt. 

In what follows, we will concentrate on two physical 
quantities: (i) the spin-squeezing parameter and (ii) the 
quantum Fisher information. 

The spin squeezing parameter, for mixed and pure 
states, is defined as m 

^2 _ N(AS’j_)min (g) 

l (^)| 2 

where N is the total atom number and (ASj_) m i n is the 
minimal variance of the spin component normal to the 

mean spin vector (S). The state is referred as the spin 
squeezed state when £ 2 < 1. 


The quantum Fisher information is an important quan¬ 
tity in interferometry. In general, the output state /5 ou t 
of an interferometer is 


Pout = e l6Sfi p in e l6Sf: 


(7) 


where Sft = S ■ ft and ft is a unit vector representing the 
effective rotation axis in a given interferometric sequence. 
The precision of the phase shift A 8 depends on the input 
state /5i n , chosen estimator and measurement performed 
on the output state [20]. According to the Cramer-Rao 
inequality, there is a lower bound on the precision with 
which the phase shift can be determined Ad >1/ y/mF Q , 
where m is the number of measurements and F Q is the 
quantum Fisher information. The quantum Fisher infor¬ 
mation, for a pure state, is given by [20] 

F q = 4(A Sl) max , (8) 

with (A^Dmax being the maximal variance of the spin 
component optimized over all possible directions ft. The 
quantum Fisher information is equal to F Q = N for the 
spin coherent state, and corresponds to the shot-noise 
limit of the phase estimation precision in optical inter¬ 
ferometry or to the projection noise in the atomic equiv¬ 
alent. The highest possible precision one can achieve is 
the Heisenberg limit A 9 > 1/yfmN, with F Q = N 2 . This 
limit can be reached using only a highly entangled state, 
for example the NOON state. 

Both quantities we are interested in are linked to the 
multiparticle entanglement. If the spin squeezing param¬ 
eter is smaller than unity, or if the quantum Fisher infor¬ 
mation is larger than the particle number, then the state 
of the system is entangled )2Q|[23l . Moreover, the quan¬ 
tum Fisher information recognizes all entangled states 
which are useful for high-sensitivity interferometry. 

Before we proceed to analyze the quantum dynamics, 
let us focus on the mean-field approximation. Although 
quantum fluctuations are lost in this description, we ar¬ 
gue that the knowledge of the classical phase space dy¬ 
namics provides an invaluable tool for the detection of 
useful states in quantum metrology. 


III. MEAN-FIELD PHASE SPACE 

The mean-field phase space dynamics is a good naviga¬ 
tor for the dynamical spin squeezing. It was shown that 
quantum evolution distinguishes between stable and un¬ 
stable classical fixed-points [21] 123, and quantum Hamil¬ 
tonian eigenstates localize on classical phase space energy 
contours ]2Sj- 

In the limit of large system size N 1, we re¬ 
place bosonic creation and annihilation operators by c- 
numbers m 

a-tVNJfre**, b-^VN^e 1 ^. (9) 

The fixed number of particles N = a 1 a + b^b, dictates 
the normalization condition p a + pb = 1. Two canonical 
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variables, the population difference z = p a — pb and the 
relative phase p = Pb ~ Pa , are sufficient to describe 
classical dynamics [2S] - Notice, when one of the modes 
is fully populated (z = 1 or z = — 1), the relative phase 
is not well defined. Spin operators become 


S x -A y\/l - Z 2 COSp, 

(10a) 

Sy -A Y \/l - z 2 sin p, 

(10b) 

N 


S z -> -j z ’ 

(10c) 


with the mean-field Hamiltonian 

N 2 /-- 

H = (H) —» — hx^p z v 1 — z2 sinp. (11) 

Equations of motion for the canonical position p and 
the conjugate momentum z can be derived from quantum 
mechanical Heisenberg equations or classical Hamilton 
equations, 



dp 

dt 

dz 

dt 


2 dll 
hNlh = 
2 dH 
HN dp 


1 — 2 z 2 

Nx vtzi? smv ’ 

(12a) 

N\z \/1 — z 2 cos p. 

(12b) 


Instead of solving these coupled differential equations 
we will analyze the topology of the phase portrait. The 
phase portrait is just a geometrical representation of tra¬ 
jectories of a dynamical system in the phase space. In our 
case, trajectories are tangent to the velocity field (p,z). 

The phase portrait in the two-dimensional system con¬ 
sists of fixed points or closed orbits. Fixed or equilibrium 
points correspond to a steady state, and satisfy 


sin p = 0 , 



(13) 


Close to a stationary point, equations of motion (121 can 
be linearized and solved exactly. Information about sta¬ 
bility of fixed points can be determined from a stabil¬ 
ity matrix (29) . Depending on eigenvalues and eigen¬ 
vectors of the stability matrix one can classify fixed 
points according to the behavior of nearby trajectories. 
The phase portrait for (11) presented in Fig. [l] consists 
of two unstable saddle fixed points located at z = 0 
and p = 0,7r, and four stable center fixed points at 
{p,z) = (±7t/2,±1/\/2). 


IV. DYNAMICS AROUND AN UNSTABLE 
FIXED POINT 

We start with the spin coherent state located at an 
unstable saddle fixed point on the equator, |7r/2,0) = 
|IV, 0) = |IV, 0) x , being the eigenstate of the S x 


FIG. I. Mean-field trajectories of the two-axis counter¬ 
twisting Hamiltonian T-L = — hx{S z S y + S y S z ). Two unstable 
saddle fixed points are located at z = 0 and p = 0, n. Four 
neutrally stable center fixed points correspond to (p, z) = 
(±7r/2, il/%/2). Colored spots visualize the spin coherent 
states around the classical unstable (red) and stable (green) 
fixed points. 


operator with eigenvalue N/2. Location of the initial 
state in corresponding mean-field phase space is sketched 
in Fig. [T] The red spot represents the initial spin coher¬ 
ent state, black lines are mean-field trajectories of the 
Hamiltonian, while the arrows indicate the direction of 
the evolution. The angle between inflowing and outflow¬ 
ing trajectories in the phase space is n/2. The regular 
dynamics takes place for short times. The initial state is 
stretched along the meridian of the Bloch sphere leading 
to the highly reduced variance of the S y component of 
the spin operator and highly increased variance of the S z 
component. The best squeezing occurs at this stage of the 
evolution, and useful states for high-precision measure¬ 
ments are generated. Moreover, the squeezing parameter 
is determined by the variance of the S y component of the 
spin operator, while the quantum Fisher information is 
determined by the variance of the S z component. Later, 
the variances of the S y and S z components of the spin 
become once again of the same order with (S x ) — —N/2. 
Next, stretching of the state takes place but the direc¬ 
tion of the evolution is opposite. The dynamics becomes 
irregular, and results depend on the total number of par¬ 
ticles. 
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A. Scaling with the system size 


In order to analyze scaling of the spin squeezing pa¬ 
rameter and the quantum Fisher information with the 
system size we use a general theory developed in mm- 
One starts with equations of motion for operators of spin 

components (Sj) which involve terms that depend on 
the first-order moment (Sj) and second-order moments 
( SjSk )■ Then, the time evolution of the second-order 
moments depends on themselves and on third-order mo¬ 
ments, and so on. It leads to the Bogoliubov-Born- 
Green-Kirkwood-Yvon (BBGKY) hierarchy of equations 
of motion for expectation values of operator products. 
We truncate the hierarchy by keeping the first- and the 
second-order moments, 

(SiSjSk) ~ (SiSjjiSk) + (SjS k ){Si) + {S k Si)(Sj) 

-2(S i )(S j )(S k ). (14) 

Let us first introduce a small parameter e = 1/N and 
transform spin components into hj = y/sSj', then the 
Hamiltonian reads 

H = * (fiyh z + khy) > (15) 

and commutation relations are [hi,hj] = ieh k Cij k . 

Equations of motion for expectation values Sj = (hj) 
and second order moments Sj k = (hjh k + h k hj) — 
2(hj)(h k ) relevant for our purposes are 

S x = 2 (Syy $zz) ? (16a) 

&yy = 4 SyySxi (15b) 

Szz = 45 ZZ S X , (16c) 


where we have introduced the dimensionless time r = 
Xt/y/e. The initial coherent state at the unstable saddle 
fixed point, |tt/ 2, 0), gives the following initial conditions: 
s x ( 0) = 1/2^/i and 6 yy (0) = <5 ZZ (0) = 1/4. Equation 
( |16a| ) takes the form s x (t) = — sinh [/ (r)] with / (r) = 
4 J 0 s x (t')cIt' , that has an analytical solution when one 
expands the function / up to the first order in Taylor 
series, /(r) ~ /(0) + /'(0)r. This requires r <C 1 or 
1/y/N. The self-consistency condition gives /(0) = 
0 and /'(0) = 4sa,(0), and the approximated solution for 
s x takes the form 


s x (t) = s x .(0) 


cosh [4s a ,(0)r] — 1 
4s x (0) 


(17) 


The above expression together with equations (16b ) and 
(16c) gives the squeezing parameter 


£2 _ g—4s x (0)r+(sinh[4s a; (0)T]— r)/4s x (0) 2 


(18) 


and the time-dependent quantum Fisher information 


F q (t) 


(0)r— (sinh[4sa, (0)r] —r)/4s x (0) 2 

£ 


(19) 


(a) 



(b) 



~ 2N 1 Y t 

>(27 tN) \ A 


FIG. 2. Scaling of the squeezing parameter (a) and the quan¬ 
tum Fisher information (b) with the total number of particles 


N. 


Notice, two approximations were made to obtain solu¬ 
tions (|l7|)-(fl9j). The first is the truncation of the BBGKY 
hierarchy, which is equivalent to the Gaussian approxi¬ 
mation, and the second is the short-time expansion which 
limits the validity of the solution to r <C 1. 

Minimization of the squeezing parameter ( |18[ ) over the 
time gives scaling of the best squeezing time Tb es t — 
In (8sa,(0) 2 ) /4s x (0) or X'tbest — ln(2N)/2N. The same 
holds for the first maximum of the quantum Fisher in¬ 
formation. Expansion of the squeezing parameter at the 
best squeezing time Tb es t in terms of the small parameter 
£ gives 



( 20 ) 


which reproduces the known result £best oc IV -1 . Leading 
terms of the maximum of the quantum Fisher informa¬ 
tion at the best time are 




^ 2 1 

best — o 5 

e £ z 


( 21 ) 


and provide the Heisenberg-like scaling F Q best — 
0.731V 2 . 

In Fig. [2] we show the squeezing parameter and the 
quantum Fisher information for different number of par¬ 
ticles N in the rescaled time. We have shown only a 
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small period of time where scaling with the number of 
particles may agree with our simulations. For further 
moments there is no time regularity. 

The scaling prediction for the squeezing parameter is 
surprisingly correct, the minima of £ 2 nicely coincide for 
different particle numbers. The prediction of the model 
for the time scaling of the first maximum of the quan¬ 
tum Fisher information is not perfect. First of all, the 
simple model predicts the position of the first maximum 
of the quantum Fisher information at the same time as 
the minimum of the squeezing parameter, since we know 
from the exact numerical calculations that the maximum 
is located at later times. Nevertheless, the time scaling 
of type xtbest ~ In (aN) / (bN) survives, and we find val¬ 
ues of a = 2n and b = 2 giving the best result. The first 
maximum of the quantum Fisher information is achieved 
at F Qib est ~ 0.677V 2 in numerical simulations, the width 
of the peak decreases with larger number of particles. 

In addition, the level of squeezing achieved by the two- 
axis counter-twisting model is always better than the 
squeezing generated by the one-axis twisting Hamilto¬ 
nian, the last scales like ~ TV -2 / 3 . The first maximum 
of the quantum Fisher information given by the two-axis 
counter-twisting model appears before the characteris¬ 
tic plateau of the quantum Fisher information generated 
by the one-axis twisting Hamiltonian which occurs at 

Xt~l/VN. 


B. Quantum states for high-precision 
measurements 


States produced by the two-axis counter-twisting 
Hamiltonian, even if they may not have a simple an¬ 
alytical form, are particularly useful for high-precision 
measurements since they provide Heisenberg-like scaling 
of the phase sensitivity A 6. Here and below, we list a 
few particular states, expanded in the Fock state basis, 
that are generated by the TACT Hamiltonian: 

(i) the Berry-Wiseman state (BW) [577] 


|BW) = 


N 

■n 


cos 


(k — N/2)n 
TV + 2 


that gives the quantum Fisher information 


I k,N-k), 

( 22 ) 


F = - 

Q 2 + TV 


N 

E 

k =0 


COS 


(k-N/2)n 

N~+~2 


(2k - TV) 2 , 


which for TV > 1 is F q « 0.137V 2 , 


(ii) the equally weighted superposition state 
(EWSS) [B] 


N 


|EWSS) = 




E \k, N ~k) 


(23) 


k=0 


with F q = (TV 2 /3) (1 + 2/TV) « 0.337V 2 , 



FIG. 3. Top panel: Time evolution of the fidelity 1261 to 
the EWSS, Twin-Fock, Berry-Wiseman and optimized over a 
Yurke states. The maxima of respective fidelity functions are: 
-Tbw, max = 0.9999 for the BW state, JEwss.max = 0.9965 
for the EWSS state, Fy, max = 0.9936 for the Yurke state 
with a = 0.678, and 7 r TF,max = 0.8732 for the twin-Fock 
state. Bottom panels: Probability distributions ( |27[ ) with 
Fock states at the maximum of the respective hdelity function 
for (a) Berry-Wiseman, (b) EWSS, (c) Yurke and (d) Twin- 
Fock states. Red crosses correspond to the ideal case given by 
equations (221-( 251. The total number of particles is TV = 50. 


(iii) the Yurke state (Y) for even number of particles [3T. 

|Y) = |TV/2 + 1, TV/2 - 1) + cos a |7V/2, TV/2) 

v 2 

+ ™|TV/2-l, TV/2 + 1), (24) 

with some real parameter a and the quan¬ 
tum Fisher information of the form Fq = 
(TV/2) ((TV/2 + 1) (2 — sin 2 a ) — 2sin 2 a), 

(iv) the twin-Fock state (TF) for even number of parti¬ 
cles (TB] 


|TF) = |TV/2, TV/2) (25) 

with F q = (TV 2 /2) (1 + 2/TV) « 0.57V 2 . 

In order to demonstrate their appearance during the evo¬ 
lution, in Fig. [3] we show fidelity functions defined as 


F A (t) = \(A\e~ mt \N,0)J 2 


(26) 
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where A is one of the states of interest (22 l-(|25|. The 


Twin-Fock and Yurke states were rotated wijth TT = 
e -i7rS x /2 £ 0 maximize j: j n addition, in Figs. |3p,-|3p we 
also show the probability distribution with Fock states 


Vk = \(k,N — k\ T) | 2 , (27) 


calculated at the maximum of respective fidelity func¬ 
tions (in gray) compared to the exact probability of the 
respective state (in red). 

The two-axis counter-twisting Hamiltonian produces 
almost ideal Berry-Wiseman, EWSS and Yurke states, 
since the maximum of the fidelity function is one, see the 
top panel in Fig. [3] Nevertheless, the probability distri¬ 
bution with Fock states is exact only for the BW state 
pTTl |32], while for the EWSS and Yurke states very small 
differences can be found. The maximum of the fidelity to 
the twin-Fock state is Ttf, max = 0.8732, showing that 
the generated state is not the perfect TF state. The 
Yurke state becomes the twin-Fock state for a = 0 what 
explains their equal fidelities T at X t = 0.065. Notice, 
the maxima of the fidelity functions for the BW, EWSS 
and Yurke states are reached at times very close to the 
best squeezing times (just before the first maximum of 
the quantum Fisher information). The maximum of the 
fidelity function for the TF state is located at times larger 
than the first maximum of the quantum Fisher informa¬ 
tion. The maxima of fidelity functions of particular states 
support the scaling of type ~ ln(aN)/(bN), where the 
coefficients a, b may depend on the state m ■ Another 
interesting observation is near the unity value of the fi¬ 
delity function for the NOON state (not shown in figures). 
Position of the maximum is not regular and shifts in time 
as the number of particles changes. 

In order to get a better insight into the regular part of 
the quantum dynamics we plot the Husimi distribution 
function 33 


Q(M,0) = \{Q,4>\ ^(t)) I 2 , (28) 

and the Wigner function W(9,(f>) j3][J at different mo¬ 
ments of time. The top panel in Fig. [4] shows a location 
of particular states on the time evolution of the quantum 
Fisher information, while the bottom panels present their 
Husimi (in the left column) and Wigner (in the right col¬ 
umn) function maps. The Husimi distribution function 
is unique, regular and positive definite. Moreover, the 
position and number of zeros of the Husimi function con¬ 
tains all information about a pure quantum state [35(536). 
This structure is visible in the logarithmic scale only un¬ 
der the quite high precision of calculations, and is invis¬ 
ible in Fig. [4] The Wigner function of a quantum state, 
in addition to localized maxima, may have interference 
fringes of negative value. We observe increasing forma¬ 
tion of the interference fringes during quantum evolution, 
even if their details may not be visible on a linear scale 
in Fig. [4] Number of fringes is equal to the total par¬ 
ticle number N. Moreover, the position of the maxima 
is the same starting from EWSS till TF states. It is 


the amplitude of the interference fringes that distinguish 
particular states. Negativity of the Wigner function was 
shown to be an indicator of non-classicality of a quantum 
state m- Thus, it is a convenient tool for detection of 
useful states for precision measurements. 

The Wigner function can be used to justify the large 
value of the quantum Fisher information. Typically, the 
presence of oriented fringes in the Wigner map increases 
its value. It will be clear when we will recover that the 
Fisher information has its geometrical interpretation |20l 
|38| f39] . If one takes 


Pout — 


—icLOSfi n. JdOSft 
Pin^ 


(29) 


with infinitesimal angle d9 then the Bures-Riemannian 
metric becomes 


^ (/4ll5 Pout) - Aq 


Pin > Sfi 


de 2 . 


(30) 


The larger the quantum Fisher information, the faster the 
state /5 out becomes distinguishable from p- m . Because the 
Wigner function is a bijective mapping we can capture 
sensitivity to rotations of a given quantum state. The 
Wigner function graphically justifies which measurement 
Sn optimizes the quantum Fisher information. If we look 
at the Wigner function of some state, it will be apparent 
that the rotation of the state through an infinitesimal 
angle dO around Sft shifts minima and maxima of inter¬ 
ference fringes making the Bures distance large. 


V. DYNAMICS AROUND A STABLE FIXED 
POINT 

The initial spin coherent state located at the classical 
stable fixed point is visualize in Fig.[l]by green spot. The 
state is slightly deformed in time leading to the Husimi 
function that changes regularly between the circular and 
elliptic shape. The frequency can be easily estimated 
within the frozen spin approximation mi mi- In this 
case it is more convenient to analyze an equivalent to |3]) 
Hamiltonian, namely, 

H = hx (s 2 x - S 2 ) , (31) 

because stable fixed points are located on the equator 

and calculations simplify significantly, H = UHW with 

JJ _ e -i(7r/2)S„ e -i(7r/4)S x _ 

In the frozen spin approximation one starts with equa¬ 
tions of motion for spin operators 

j t S x =x(s y S z + S z Sy), (32a) 

^S v = -x(sJ x +Sj x ), (32b) 

j t S z = 2 X ( S X Sy + S V S X ) . (32c) 
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FIG. 4. Top panel: arrows indicate the locations of quantum states that are plotted in bottom panels on the time evolution of 
the quantum Fisher information. Bottom panels: the Husimi (left column) and Wigner (right column) functions at particular 
moments of time: (A) the initial coherent spin state, (B) the maximum of the BW fidelity function, (C) the maximum of the 
EWSS fidelity function, (D) the best spin squeezing, (E) the maximum of the Yurke fidelity function, (F) the maximum of the 
quantum Fisher information, (G) the maximum of the twin-Fock fidelity function, and (H) the second local minimum of the 
quantum Fisher information. The total number of particles is N = 50. 















































(a) 



\f2~Nxt 


FIG. 5. The spin squeezing parameter (a) and the quantum 
Fisher information (b) for different numbers of particles N 
around a stable fixed point. 


The initial condition |7r/2,0) corresponds to ( S x ) = N/2 
and the dynamics will be captured around this point. In 
the frozen spin approximation one replaces the operator 
S x by its mean value ( S x ) = N/2, and ends up with 
equations for the remaining operators 

J t Sy = -xNS z , (33a) 

j t S z = 2 X NS V . (33b) 

One can easily solve these coupled differential equations 


the optimized quantum Fisher information are 

= "i-isinV), (35a) 

Fq = 4(AS 2 ) = N [1 + sin 2 (wf)] . (35b) 

Results of the frozen spin approximation can be com¬ 
pared to the exact numerical simulations presented in 
Fig. [5] Although overall shape of the function slightly 
deviates from sin 2 (wt), extreme values and the frequency 
are in a good agreement. Now it is clear why the ini¬ 
tial spin coherent states located at classical stable fixed 
points give shot-noise limited sensitivity for quantum 
metrology. The maximal value of the Fisher informa¬ 
tion scale linearly with the particle number N. It is not 
possible to beat shot-noise limit in this configuration. 

VI. SUMMARY 

We have discussed the quantum dynamics generated 
by the two-axis counter-twisting Hamiltonian in the con¬ 
text of quantum metrology. We have started the anal¬ 
ysis within the mean-field description which identifies 
a convenient location of the initial spin coherent state. 
The quantum dynamics around a stable fixed point sup¬ 
ports the shot-noise-limited sensitivity for precision mea¬ 
surements. On the other hand, the quantum dynam¬ 
ics around an unstable fixed point generates quantum 
states that give Heisenberg-like scaling. We have cal¬ 
culated the spin squeezing parameter and the quantum 
Fisher information explaining their scaling with the sys¬ 
tem size. Our results show that the quantum dynamics 
with the two-axis counter-twisting Hamiltonian creates 
quantum correlations in a regular way on a short time 
scale. In addition, a characteristic feature of the output 
of TACT Hamiltonian are concentric fringes of negative 
values in the Wigner functions of states located around 
the first maximum of the quantum Fisher information. 
It may be difficult to identify an efficient interferomet¬ 
ric strategy and measurement to fully recover this part 
of the quantum Fisher information. In the OAT model, 
for instance, the quantum Fisher information can be re¬ 
trieved by a standard measurement of the population im¬ 
balance |32] 05]. 


S y {t) = S y (0) cos(wf)-^^-(O) sin(wt), (34a) 

v 2 

S z (t) = S'z(O) cos(wt) + V2S y (0) sin(wt), (34b) 
with w = \[2N\. The spin squeezing parameter £ 2 and 
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